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ABSTRACT 



Context. Giant planets orbiting main-sequence stars closer than 0.1 AU are called hot Jupiters. They interact with their stars affecting their 
angular momentum. 

Aims. Recent observations provide suggestive evidence of excess angular momentum in stars with hot Jupiters in comparison to stars with 
distant and less massive planets. This has been attributed to tidal interaction, but needs to be investigated in more detail considering also other 
possible explanations because in several cases the tidal synchronization time scales are much longer than the ages of the stars. 
Methods. We select stars harbouring transiting hot Jupiters to study their rotation and find that those with an effective temperature ^ 6000 K 
and a rotation period P lot 5 10 days are synchronized with the orbital motion of their planets or have a rotation period approximately twice that 
of the planetary orbital period. Stars with T c k S 6000 K or P m[ ^ 10 days show a general trend toward synchronization with increasing effective 
temperature or decreasing orbital period. We propose a model for the angular momentum evolution of stars with hot Jupiters to interpret 
these observations. It is based on the hypothesis that a close-in giant planet affects the coronal field of its host star leading to a topology with 
predominantly closed field lines. An analytic linear force-free model is adopted to compute the radial extension of the corona and its angular 
momentum loss rate. The corona is more tightly confined in F-type stars, and in G- and K-type stars with a rotation period shorter than ~ 10 
days. The angular momentum loss is produced by coronal eruptions similar to solar coronal mass ejections. 

Results. The model predicts that F-type stars with hot Jupiters, r c n ^ 6000 K and an initial rotation period S 10 days suffer no or very little 
angular momentum loss during their main-sequence lifetime. This can explain their rotation as a remnant of their pre-main-sequence evolution. 
On the other hand, F-type stars with P TOt > 10 days, and G- and K-type stars experience a significant angular momentum loss during their 
main-sequence lifetime, but at a generally slower pace than similar stars without close-in massive planets. Considering a spread in their ages, 
this can explain the observed rotation period distribution of planet-harbouring stars. 

Conclusions. Our model can be tested observationally and has relevant consequences for the relationship between stellar rotation and close-in 
giant planets as well as for the application of gyrochronology to estimate the age of planet-hosting stars. 

Key words, stars: planetary systems - stars: late-type - stars: rotation - stars: magnetic fields - stars: coronae 
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The search for planetary systems has revealed a population of 
planets having a mass comparable to that of Jupiter orbiting 
main-sequence late-type stars closer than ~ 0.1 AIjQ. They in- 
clude about 25 percent of all known planets and have been 
called hot Jupiters because of their close proximity to their host 
stars. Remarkable interactions between those planets and their 
host stars are expected, both as a consequence of tides or recon- 
nection between planetary and stellar magnetic fields (cf., e.g., 
Cuntz et al 2000; Gu & Suz uki 2009). Eviden c e of a magnetic 



quences for the evolution of the stellar angular momentum. 
The presence of a close-in giant planet may significantly af- 
fect the structure and the energy balance of the coronal field. 
Kas hvap et al.1 (120081) show that stars with hot Jupiters have X- 
ray luminosities up to 3 - 4 times greater than stars with distant 
planets. This suggests that a close-in planet may enhance mag- 
netic energy dissipation or lead to a predominance of closed 
and brighter magnetic structures in a stellar corona. 



interaction was first shown by IShkolnik et all d2003l) and ha s 



LanzaN2008l) . lLanzal (§009) 



bee n recently reviewed b y, e.g 
and lShkolnik et al.1 d2009l) . 

In the present study, we focus on the modifications induced 
in a stellar corona by a close-in giant planet and their conse- 
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Pont ( 20091) has found that the sample of stars with tran- 
siting hot Jupiters show a statistic excess of rapidly rotat- 
ing objects in comparison to stars without close-in planets. 
Plotting stellar rotation rate vs. the orbital semimajor axis, 
normalized to the average of the stellar and planetary radii, 
and the planet-to-star mass ratio, he finds some empirical ev- 
idence of faster rotation in stars with closer and more mas- 
sive planets. This may be interpreted as an indication of tidal 
interaction that drives stellar rotation toward synchronization 
with the planetary orbital period because the tidal torque is ex- 
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pected to increase for closer and more massive planets (e.g., 
Mardling & L~inll2002l) . 

There are a few examples of stars whose rotation appears 
to be synchroni zed with the orbit of their close-in planets, no- 
tably t B ootis (Donatietal. 2008|). and the tran siting system 



CoRoT-4 (lAigrain et al.ll2008t lLanza et al.ll2009h . They are F- 
type stars with a shallow outer convective envelope. In par- 
ticular, t Boo has an estimated age of « 2 Gyr and is orbited 
by a planet with Msin; = 4.38 My corresponding to » 6.8 
Jupiter masses at a distance of ~ 0.049 AU, i.e., ~ 7.2 stellar 
radii, if an inclination i 40° is adopted, as suggested by stel 



lar D oppler Imaging models (ILeigh et al. I 12003c ICatala et al 
2007). Assuming that only the envelope of t Boo is in a syn- 



chronous rotation state. iDonati et al. (2008) found a synchro- 
nization timescale compatible with the main-sequence lifetime 
of the star. However, in the case of CoRoT-4a that timescale 
is of the order of 350 Gyr, because the semimajor axis of 
the planetary orbit is 17.4 stellar radii and the mass of the 
planet is only 0.72 Ju piter masses mak ing any tidal interac- 
tion extremely small (lLanza et al.l 120091) . Therefore, a differ- 
ent process is required to account for the synchronization of 
CoRoT-4a. Two other intriguing cases are those of the host 
stars of XO-4 and HAT-P-6 whose rotation periods are close 
to twi ce the orbital periods of their transiting planets, respec- 
tively dMcCullough et al]|2008l) . 

Tidal interactions are not the only processes affecting spin 
and orbital angular momenta. Since late-type stars have mag- 
netized stellar winds that produce a remarkable braking of their 
rotation during their main-sequence lifetime, a continuous loss 
of angular momentum must be taken into account to model 
the evol ution of the stellar s pin in a proper way. As shown 
by, e.g., iDobbs-Dixon et al.l (120041) . this may have a signifi- 



cant impact not only on the evolution of stellar rotation but 
also on that of the orbital parameters, notably the eccentricity 
and the semimajor axis, especially during the initial stages of 
the main-sequence evolution when the star is a fast rotator and 
stores most of the angular momentum of the system that can 
be transferred to the planetary orbit to excite its eccentricity 
an d/or inc r ease i ts semimajor axis. 

Lanzal ( 200 8) proposed a model for the interaction between 



the coronal and the planetary magnetic fields considering a lin- 
ear force-free equilibrium for the coronal field. Here we shall 
apply that model to study the angular momentum loss from the 
coronae of stars hosting hot Jupiters assuming that a close-in 
planet leads to a corona with predominantly closed magnetic 
field lines because it tends to reduce the magnetic helicity of 
the field via a steady dissipation of magneti c energy a ssociated 
with its motion through the corona (cf. lLanzal 120091) . Fields 
with a lower helicity are characterized by a topology with a 
greater fraction of closed field lines, while an increase of the 
helicity beyond a certain threshold may lead to an openi ng up 
of the field lines dZhang et al.ll2006l;IZhang & Flverj|2008l) . The 
recent MHD simulations of the effects of a close-in massive 



planet on the coronal field of a star by Cohen et a l. (2009) con- 
firm that the planet inhibits the expansion of the coronal field 
and the acceleration of the stellar wind. 

Under the hypothesis of a coronal field with predominantly 
closed field lines, we find that a close-in planet may signifi- 



cantly reduce magnetic braking in rapidly rotating stars. This, 
in addition to tidal interaction, may explai n the tendency to- 
ward synchronization found by IPond (12009b if the initial ro- 
tation periods of the stars are close to the orbital periods of 
their planets, in agreement with models of rotational evolution 
which assume that both a star and its planet are dynamically 
coupled to a circumstellar disc during the first few million years 
of their evolution (cf. Sect. 14. 11 1. 



2. Properties of transiting planetary systems 

The rotation periods of planet-hosting stars are difficult to mea- 
sure from the rotational modulation of their flux because plan- 
ets can be detected by radial velocity monitoring only around 
relatively inactive stars, i.e., with a photometric modulation be- 
low sa 0.01 mag. The situation is going to change with space- 
borne phojxjmej^maj^ and Kepler (cf.. 
e.g.. " ' """" 



Alonso et al. 2008; Aigrain et al. 2008), but at the mo 



ment the most reliable measurements of the rotation rates of 
planet-harbouring stars come from spectroscopy, viz., from ro- 
tational line broadening v sin i. To find the rotation velocity v, 
we need to know the inclination of the rotation axis i, which 
can be derived in the case of transiting systems by assuming 
that the stellar and orbital angular momenta are aligned. This 
alignment results from the formation of planets in a circumstel- 
lar disc, if the gravitational interaction among planets is not too 
strong, otherwise one expects a significant misalignment of the 
spin and orbit al angular momenta a s well as eccentric planetary 
orbits (e.g., iNagasawa et al . 2008). It is interesting to note that 
the angle A between the projections of the spin and orbital an- 
gular momenta on the plane of the sky can be measured through 
the Rossiter-McLa ughlin effect that has already been detecte d 
in several systems dOhta et alj|2005l iFabrvckv & Winnll2009l) . 

In view of the possibility of guessing the value of the in- 
clination, we focus on transiting systems. As to the end of 
September 2009, 62 transiting systems are knowr@. We es- 
clude those with an eccentricity e > 0.06 and/or with a sig- 
nificant misalignment as indicated by the Rossiter-McLaughlin 
effect. They are listed in Table Q] where we report, from the 
first to the fourth column, the name of the system, its or- 
bital eccentricity, the misalignment angle A as derived from 
the Rossiter-McLaughlin effect, and the references, respec- 
tively. Moreover, we exclude CoRoT-7 bec ause its planets are 
not hot Jupiters, bu t terrestrial-sized objects dLeger et al.l l2009; 
Oueloz et al.ll2009l) : CoRoT-6 because it was announced, but 
its parameters have not been published yet; WASP-2, WASP-8, 
OGLE-TR-111, OGLE-TR-182, and OGLE-TR-211 because 
no data on their stellar rotation were found in the literature. 

The relevant parameters of the considered systems are 
listed in Tables [2] and [3] In Table [2] from the first to the eighth 
column, we list the system name, the orbital period P 01 -b, the 
semimajor axis of the planetary orbit a, the radius of the planet 
R p , the mass of the planet M p , the radius of the star R, the mass 
of the star M, and its effective temperature 7" e ff, respectively; in 
Table [3] from the first to the seventh column, we list the sys- 
tem name, the v sin i of the star, the rotation period of the star 
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Fig. 1. Upper panel: The synchronization parameter n/Q vs. the effective temperature of the star in transiting planetary systems; 
lower panel: an enlargement of the lower portion of the upper panel, to better show the correlation close to n/Q = 1 and n/Q = 2, 
which are marked by horizontally dashed lines. The names of the systems are reported in both panels, although they are omitted 
for n/Q < 7.5 in the upper panel to avoid confusion. 



as derived from the rotational modulation of its flux P rot , the 
minimum and maximum estimated age for the star, the syn- 
chronization time for stellar rotation r sync , the eccentricity of 
the planetary orbit e, and the references for both Tables [2] and 
[3] respectively. When it is not determined from the rotational 
modulation of the flux, we estimate the rotation period of a star 
as: P TOt = 2nR/(v sini), where the radius of the star is taken 
from Table [2] and an inclination i = 90° is adopted, which is 
perfectly adequate for transiting systems in view of the aver- 
age precision of our v sin i and stellar radii. In this case, we 



do not list F rot i n Tab le [3] We make use of the formulae of 
Mardling & Lin (2002) to compute r sync as: 



"*~sync 



l 

Q 



da 



dt 



1 



2 y 2 Q' t \ M 



Mr, 



(§r-(s) 



GM 



(1) 



where yR =t 0.35/? is the gyration radius of the star, <2» is 
the modified tidal quality factor of the star defined as Q' t = 
3Q m /2k t , where g* is the specifi c dissipation function an d 



k* is the Love number (cf., e.g., iMurrav & Dermott Il999h : 
n = 27r/P or b is the mean orbital motion, 



Q = 2jr/P ml the an- 
gular velocity of rotation of the star, and G the gravitation con- 
stant. Eq. ([TJ is valid for circular orbits and when the spin axis 
is aligned with the orbital angular momentum. To compute the 
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Fig. 2. The synchronization parameter n/Q vs. the stellar rotation period f rot for our sample of transiting extrasolar planetary 
systems. Different symbols indicate different ranges of stellar effective temperatures: filled circles: r e ff > 6500 K; open circles: 
6000 < r e jf < 6500 K; triangles: 5500 < T e ff < 6000 K. The vertically dotted line marks the rotation period of 10 days below 
which systems show a tendency to cluster around the values n/Q — 1 (synchronization) and n/Q = 2. 



Name 



References 



HD 17156 
HD 80606 

HD 147506/HAT-P-2 

GJ436 

XO-3 

HAT-P-7 

HAT-P-11 

WASP- 14 

WASP- 17 

CoRoT-1 



0.6819 + 0.0044 



Barbieri et al. (2009 



0.9332 ± 
0.5163 - 
0.160 ± 
0.2884 ± 

0.198 ± 
0.091 ± 
0.129 



0.0008 
: 0.025 
0.019 
0.0035 

0.046 
0.003 
±0.1 



Fischer & Valenti (2005); Pont et al. (2009a) 
Bakos et al. (2007a) : Loeillet et al. (2008) 
Bean et al. (2008) : Demorv et al. (2007) : Maness et al. 
37° ± 4° Johns-Krull et al. (2008) : Winn et al. (2009a) 
182° ±10° Narita et al. (2009) : Winn et al. (2009b) 

Dittmann et al. (2009) : Bakos et al. (2009c) 
-33° ±8° Joshi et al. (2009) : Johnson et al. (2009b) 
-147° ±50° Anderson et al. (2009) 
77° ±11° Pont et al. (2009a) 



Table 1. Transiting planetary systems excluded from the present analysis. 



synchronization times listed in Table|3] we adopt Q' t = 10 6 and 
assume that the entire star is synchronized as it is customary 
in tidal theory and is suggested by the tidal evolution of close 
binaries observed in stellar clusters of different ages. If n = Q, 
i.e., when the stellar rotation is synchronized with the orbital 
motion, there is no tidal contribution to the variation of the stel- 
lar angular momentum. If we consider the case of F-type stars 
discussed in the Introduction, i.e., t Boo and CoRoT-4a, since 
their initial rotation is synchronized with the orbital motion, 
tidal effects can be neglected. Nevertheless, we list in Table [3] 
a lower limit for the synchronization timescale for CoRoT-4a 
computed by assuming Q = 2n. 



We note that r sync regards the angular momentum exchange 
between the orbital motion and the stellar spin. On the other 
hand, orbital circularization proceeds much faster owing to the 
dissipation of energy inside the planet at almost constant orbital 
angular momentum. The dissipation of energy inside the star is 
smaller by a factor of 10~ 2 - 10~ 3 when we adopt Q' v — 10 s for 
the planet and Q'„ = 10 6 (cf., e.g., iMatsumura et all2008l) . 

As a measure of synchronization of a planet-harbouring 
star, we adopt the ratio n/Q = P ro t/P bs- Its dependence on 
the orbital p ara meters has already be en investigated by, e.g., 
IPonj d2009l) and iLevrard et al.1 d2009l) . Here, we focus on the 
correlation between n/Q and the effective temperature of the 
star that is plotted in Fig.Q] Although there is a large scatter for 
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T e ff ^ 6000 K, there is a general trend toward synchronization, 
i.e., a decrease of n/Q toward unity, with increasing effective 
temperature and the synchronized systems have T e ff > 6200 K. 
We note that for 6000 < T ef( < 6500 K there is still a signif- 
icant scatter in n/Q. Nevertheless, two subgroups of systems 
can be identified: one consists of those being close to n/Q = 1 
or to n/Q = 2, while the other consists of those showing n/Q 
remarkably greater than 2. These two groups appear to be sep- 
arated if we plot n/Q vs. P rot , as it is shown in Fig.|2j with all 
systems with P lot < 10 days having n/Q 2, with the excep- 
tion of WASP- 18. 

The significance of the clustering of the systems around 
n/Q = 1, 2 can be assessed by means of the Kolmogorov- 
Smirnov test (hereinafter KS) based on the cumula tive distri- 
bution function of the observed values (see 



e.g., 



Press et al 



19921) . In the upper panel of Fig. [3] we compare the cumulative 
distribution functions for the systems with T e ff > 6000 K and 
Teff < 6000 K with the corresponding uniform distributions for 
n/Q < 12, respectively. The probability that the observed dis- 
tribution of n/Q for r e ff > 6000 K is drawn from a uniform 
distribution is only 2.50 x 10~ 6 , according to the KS statistics. 
On the other hand, the probability that the distribution of n/Q 
for r e ff < 6000 K is compatible with a uniform distribution is 
0.575. In the lower panel of Fig. [3] we plot an histogram with 
bin sizes of unity to better show the clustering of n/Q around 
the values < 2 for T e g > 6000 K and the almost uniform dis- 
tribution for r e jf < 6000 K. In Table |4] we list the results of 
the KS test for different maximum values of n/Q to show that 
the clustering around the values 1 and 2 for > 6000 K 
is significant and does not depend on the choice of the upper 
limit of n/Q used to define the data sample. Specifically, in the 
columns from the first to the fifth, we list the maximum value 
of n/Q defining the sample, the probability Pi that the subset 
with T e ff > 6000 K is drawn from a uniform distribution, the 
number N\ of systems in that subset, the probability P2 that the 
subset with r e jf < 6000 K is drawn from a uniform distribution, 
and the number No of systems in that subset, respectively. For 
subsets containing less than 6 systems, we do not give the cor- 
responding probability because the KS statistics is computed 
by a n asymptotic fo rmula which is not accurate in these cases 
(cf. IPress et al.|[l992l) . 

The trend toward synchronization seen in Fig. Q] could be 
due to the general decrease of the rotation periods of main- 
sequence stars with increasing effective temperature and the 
clustering of the orbital periods of the hot Jupiters around 
~ 3 - 4 days, irrespective of any kind of star-planet interac- 
tion. To test this explanation, we compare the evolution of the 
rotation periods of planet-hosting sta rs with those o f the stars 
without planets as parameterized by iBarnesI (12007b from the 
study of stellar rotation in open clusters and in the field. It is 
important to consider the stellar age in addition to the effec- 
tive temperature because stellar rotation periods on the main 
sequence show a re markabl e depen dence on both those param- 
eters. According to Barnes d2007h . the rotation period P mt in 
days and the age t in Myr are related to the B - V colour index 
according to the formula: 



where n = 0.5189 + 0.007 is almost identical to the value 1/2 
of the well-known Skumanich braking law, a = 0.7 725 + 0.01 1 
and b = 0.601 + 0.024 (cf. Eq. (3) in lBarneslE007l) . It is inter- 
esting to note that Eq. with n = 1/2, i.e., the Skumanich's 
law, can be obtained by assuming an angular momentum loss 
rate: dL/dt = K(B-V) Q 3 , where K oc a[(B-V)-0A] b depends 
on the parameters a and b of the Barnes' relationship. 

In view of the dependence on the stellar age, we restrict our 
comparison to the 24 systems with an age estimate in the liter- 
ature and plot in Fig.|4]their P ro t*~" vs - T e s, where n = 0.5189 
and the effective temperature has been con verted into the B — V 
colour index thanks to the calibration by BesselJ dl979l) . The 



P mt r" =a[(B-V)-0A] b 



(2) 



solid line is the correlation found by Barnes (2007) for stars 
without hot Jupiters. The two systems with P e ff > 6700 K are 
CoRoT-3 and OGLE-TR-L9 whose rotation periods are signif- 
icantly affected by tidal interaction (see below), so they are ex- 
cluded from our analysis. 

Ax 2 test of the goodness of fit of the Barnes' relationship 
for the remaining 22 systems has been performed and gives a 
probability of 0.14 that the obtained^ 2 = 28.96 is compatible 
with Eq. ©, suggesting that stars with transiting hot Jupiters 
are, on the average, faster rotators than similar stars without 
planets of the same age. The goodness-of-fit probability de- 
creases to 0.018 if we restrict the comparison to the 15 stars 
with 5500 < r e ff < 6700 K Of 2 = 28.61), reinforcing the con- 
clusion for this subsample of stars. 

On the other hand, assuming a = 0.56, we obtain 
X 1 = 12.66 for the 22 systems with P e ff < 6700 K, which 
has a goodness-of-fit probability of 0.94, indicating that a 
Skumanich-type law with a reduced angular momentum loss 
rate is an adequate description of the overall evolution of the 
rotation of stars with hot Jupiters. This implies that the rota- 
tion periods of the planet-hosting stars are, on the average, a 
factor of 0.7 shorter than those of the stars without planets of 
the same age. This is particularly significant because the detec- 
tion of transits and the radial velocity confirmation of planets 
introduce a bias toward stars with a lower level of activity and 
therefore longer rotation periods. An important consequence of 
the lower angular momentum loss rate found in stars with hot 
Jupiters is that the tendency toward synchronization shown in 
Fig.[T]cannot be explained on the basis of the general decrease 
of the rotation period with effective temperature observed in 
main-sequence stars without hot Jupiters. 

The general trend shown in Fig.[T]as well as the correlation 
in Fig. 12 may result from tidal effects, but a plot of the ratio 
Tsync/ T evoi> where T evo i is an estimated upper bound for the stel- 
lar age, vs. the effective temperature, shows a general increase 
of that ratio with the effective temperature (cf. Fig. [5]). Morover, 
~ 70 percent of the systems for which we have an age esti- 
mate have not had enough time to synchronize the rotation of 
their stars. The value of r sync is one or two orders of magnitude 
greater than the maximum estimated stellar age in the case of 
CoRoT-4, HAT-P-9 (n/Q * 1), HAT-P-6, or XO-4 (n/Q 2). 
Even by decreasing the value of Q' t by one order of magnitude, 
we cannot eliminate the discrepancy. Therefore, tidal effects 
alone do not appear to be a viable explanation for the synchro- 
nization observed in several systems. A similar conclusion is 
reached for the correlation between n/Q and P lot seen in Fig. [2] 
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Fig. 3. Upper panel: Cumulative distribution functions C of 
n/Q for the systems with a stellar effective temperature T e g > 
6000 K (black solid line) and r eff < 6000 K (red dashed line) 
and n/Q < 12; the expected cumulative distributions for a uni- 
form distribution of n/Q are also plotted for T e ff > 6000 K 
(black dotted line) and T e g < 6000 K (red dot-dashed line), 
respectively. Lower panel: Histogram of the distribution of 
n/Q for systems with > 6000 K (black solid line) and 
T e g < 6000 K (red dashed line), respectively. 




Fig. 4. Age-normalized rotation period, P lot r" with n = 
0.5189, vs. the stellar effective temperature T e g for our sam- 
ple ofstos with transiting hot Jupiters. The relationship found 
by iBarnesI (12007b for stars without hot Jupiters is plotted with 
a solid line. The dashed line plots the same relationship pro- 
posed by Barnes but for a lower angular momentum loss rate 
parameter a = 0.56. 
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Table 4. Kolmogorov-Smirnov test of uniform distribution for 
different subsamples of the distribution of n/Q.. 



because systems such as CoRoT-4, HAT-P-6, HAT-P-9, orXO- 
4 have too long synchronization time scales to explain their 
low values of n/Q.. This implies that another mechanism must 
be at work, in addition to tidal effects, to produce the observed 
trend toward synchronization with increasing effective temper- 
ature and the remarkable concentration of systems around the 
values n/Q = 1 and n/Q - 2 observed for P rot < 10 days 
and r e ff > 6000 K. We conjecture that such a mechanism is 
related to the modification of the stellar coronal field and to 
the different magnetically-controlled angular momentum loss 
induced by a close-in massive planet. In the framework of such 
a hypothesis, we introduce in the next Sections a model for the 
angular momentum content of the corona of a star harbouring 
a hot Jupiter and discuss how the evolution of its angular mo- 
mentum is affected. 



3. Modelling the angular momentum content of a 
stellar corona 

3.1. Coronal field model 

We adopt a spherical polar coordinate frame having its origin 
at the barycentre of the host star and the polar axis along the 
stellar rotation axis. The radial distance from the origin is indi- 
cated with r, the colatitude measured from the North pole with 
8, and the azimuthal angle with <p. The planet orbit is assumed 
circular and lying in the equatorial plane of the star according 
to the selection criteria considered in Sect. [2] We adopt a refer- 
ence frame rotating with the angular velocity of the star Q with 
respect to an inertial frame. 

We model the coronal field under the hypothesis that the 
magnetic pressure is much greater than the plasma pressure and 
the gravitational force, so we can assume a force-free magne- 
tohydrostatic balance, i.e., the current density J is everywhere 
parallel to the magnetic field B, viz. J x B = 0. This means 
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Fig. 5. The ratio between the tidal synchronization timescale r sync and an estimated upper limit for the age of the star T evo i vs. the 
effective temperature for the systems for which an estimate of the stellar age has been found in the literature (cf. Table|3]l. 



that V x B — aB, with the forc e-free parameter a constant 
along each field line (|Priesllll982i) . If a is uniform in the stellar 
corona, the field is called a linear force-free field and it satisfies 
the vector Helmoltz equation V 2 B + a 2 B = 0. Its solutions in 
spheri cal g eometry have been studied by, e. g 
( 1956) and lChandrasekhar & Kendalil dl957l) . 



Chandrasekhar 



Linear force-free fields are particularly attractive in view of 
their mathematical symplicity and their minimum-energy prop- 
erties in a finite domain that contains all the lines of force, as 



shown by Woltier ( 1958). Specifically, in ideal magnetohydro- 
dynamics, the minimum energy state of a magnetic field in a 
finite domain is a linear force-free state set according to the 
boun dary conditi ons and the conservation of magnetic helic- 
ity. iBergerl (11985b modified the definition of magnetic helicity 
introducing a relative magnetic helicity which is conserved in 
spite of the fact that the lines of force may not be contained 
into a finite volume. Such a relative helicity is the relevant con- 
served quantity in the case of a stellar corona whose lines of 
force cross the surface of the star. 

Considering th e cons ervat ion of relative helicity , e.g., 
Hevvaerts & Priestl d 1984b and iRegnier & Priest! d2007l) have 
conjectured that the minimum energy state actually allowable 
to a solar active region in a semi-infinite atmo sphere is a lin- 
ear fo rce-free state (see, however, discussion in IZhang & Low 
20051) . We conjecture that the enhanced dissipation induced by 
the motion of a hot Jupiter inside a stellar corona drive s the 
coron al field toward such a linear force-free state (cf. Lanzal 
2009). Note that, even in the case of a synchronous system 



such as t Boo or CoRoT-4, the differential rotation of the stel- 
lar surface implies relative velocities of the order of 10 - 30 km 



s between the planet a nd the coronal field (cf., Ca tala et al 



2007; Lan za et al.l l2009). 



To model the stellar coronal field, we consider only the 
dipole-like component (i.e., with a radial order k — 1) of the 
linear force-free solution of IChandrasekhar & Kendalil (11957b 
because it has the slowest decay with distance from the star 
and therefore contains most of the angular momentum of the 
corona (see below and Sect. [33}. Moreover, an axisymmetric 
field (i.e., with an azimuthal degree m = 0) is the simplest ge- 
ometry to model the corona and is also favoured by models of 



magnetic star-planet interaction, as discussed by Lanza (2008, 
2009b . In particular, an axisymmetric field reproduces the ob- 



served phase lag between the planet and the chromospheric hot 
spot i nduced by its interaction wi th the coronal field of the star 
(e.g., IShkolnik et ai1l2005[l2008l) . 

Our li near force- f ree fie ld can be expressed in the for- 
mulism of Fiver et al.l (120041) as: 



B 



1 



r sin 6 



\dA dA. 
~r~dB r ~~d~r 



0+Q(A)4, 



(3) 



where A(r, 6) is the flux function of the field and Q = aA. 
Magnetic field lines lie over surfaces of constant A(r, 6), as 
follows by noting that B ■ VA = 0. The flux function for our 
dipole-like field geometry is A(r, 0) = B(,R 2 g(q) sin 2 6, where 
2Bo is the magnetic field intensity at the North pole of the star, 
R the star's radius and the function g(q) is defined by: 



8(9) = 



[b J-3f 2 (q) + coh/ 2 (q)] yfq 
[b J-3 /2 (qo) + c J3/2(qo)] y/qo' 



(4) 
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Name P orb a R p M p R M T cff 

(days) (AU) (R s ) (Mj) (i? G ) (M Q ) (K) 



HD 149026 


2 


.87589 


0.04320 


0.728 


+ 


0.012 


0.360 ± 


0.030 


1.497 


± 


0.069 


1.300 


± 


0.060 


6147 ± 


:50 


HD 189733 


2 


,21858 


0.03100 


1.154 


+ 


0.017 


1.150 ± 


0.040 


0.755 


+ 


0.011 


0.820 


+ 


0.030 


5050 + 


: 50 


HD197286/WASP7 


4 


,95466 


0.06180 


0.925 


+ 


0.043 


0.960 ± 


0.200 


1.236 


± 


0.050 


1.280 


+ 


0.160 


6400 ± 


100 


HD209458 


3 


,52475 


0.04700 


1.320 


± 


0.025 


0.657 ± 


0.006 


1.125 


± 


0.022 


1.101 


+ 


0.064 


6117 ± 


:26 


TrES-1 


3 


,03007 


0.03930 


1.081 


+ 


0.029 


0.760 ± 


0.050 


0.811 


+ 


0.020 


0.890 


± 


0.035 


5250 + 


200 


TrES-2 


2 


,47063 


0.03670 


1.221 


+ 


0.044 


1.198 ± 


0.053 


1.000 


± 


0.035 


0.980 


± 


0.062 


5850 + 


: 50 


TrES-3 


1 


,30619 


0.02280 


1.336 


+ 


0.034 


1.910 ± 


0.080 


0.830 


+ 


0.020 


0.930 


+ 


0.040 


5650 + 


: 75 


TrES-4 


3 


,55395 


0.05100 


1.783 


+ 


0.090 


0.925 ± 


0.080 


1.849 


± 


0.090 


1.400 


± 


0.100 


6200 + 


: 75 


XO-1 


3 


,94153 


0.04880 


1.189 


± 


0.023 


0.900 ± 


0.070 


0.928 


+ 


0.033 


1.000 


± 


0.030 


5750 ± 


: 13 


XO-2 


2 


,61584 


0.03680 


0.984 


± 


0.019 


0.570 ± 


0.060 


0.964 


+ 


0.020 


0.980 


+ 


0.030 


5340 + 


: 32 


XO-4 


4 


,12502 


0.05550 


1.340 


± 


0.050 


1.720 ± 


0.200 


1.560 


± 


0.050 


1.320 


+ 


0.020 


6397 + 


: 70 


XO-5 


4 


,18775 


0.04880 


1.089 


+ 


0.057 


1.077 ± 


0.037 


1.080 


+ 


0.050 


0.880 


+ 


0.030 


5370 + 


: 70 


HAT-P-1 


4 


,46529 


0.05510 


1.203 


+ 


0.051 


0.530 ± 


0.040 


1.115 


± 


0.043 


1.120 


+ 


0.090 


5975 + 


:45 


HAT-P-3 


2 


,89970 


0.03894 


0.890 


+ 


0.046 


0.599 ± 


0.026 


0.824 


+ 


0.040 


0.936 


+ 


0.050 


5185 + 


:46 


HAT-P-4 


3 


,05654 


0.04460 


1.270 


± 


0.050 


0.680 ± 


0.040 


1.590 


+ 


0.070 


1.260 


± 


0.100 


5860 + 


: 80 


HAT-P-5 


2 


,78849 


0.04075 


1.260 


± 


0.050 


1.060 ± 


0.110 


1.170 


+ 


0.050 


1.160 


± 


0.060 


5960 ± 


100 


HAT-P-6 


3 


,85298 


0.05235 


1.330 


± 


0.061 


1.057 ± 


0.119 


1.460 


+ 


0.060 


1.290 


± 


0.060 


6570 + 


: 80 


HAT-P-8 


3 


,07632 


0.04870 


1.460 


+ 


0.020 


1.520 ± 


0.170 


1.580 


± 


0.070 


1.280 


+ 


0.040 


6200 + 


: 80 


HAT-P-9 


3 


,92289 


0.05300 


1.400 


+ 


0.060 


0.780 ± 


0.090 


1.320 


+ 


0.070 


1.280 


+ 


0.130 


6350 + 


150 


HAT-P-10/WASP11 


3 


,72247 


0.04350 


1.055 


+ 


0.040 


0.487 ± 


0.030 


0.790 


± 


0.030 


0.830 


+ 


0.030 


4980 + 


: 60 


HAT-P-1 2 


3 


,21306 


0.03840 


0.963 


± 


0.025 


0.211 ± 


0.012 


0.700 


+ 


0.030 


0.730 


± 


0.020 


4650 + 


: 60 


HAT-P-1 3 


2 


,91626 


0.04260 


1.280 


± 


0.080 


0.851 ± 


0.046 


1.560 


+ 


0.080 


1.220 


± 


0.100 


5638 + 


: 90 


WASP-1 


2 


,51996 


0.03820 


1.443 


± 


0.039 


0.790 ± 


0.130 


1.453 


± 


0.032 


1.150 


± 


0.090 


6110 ± 


245 


WASP-2 


2 


,15223 


0.03070 


1.038 


+ 


0.050 


0.880 ± 


0.070 


0.780 


+ 


0.060 


0.790 


+ 


0.150 


5200 ± 


200 


WASP-3 


1 


,84684 


0.03170 


1.255 


+ 


0.085 


1.760 ± 


0.090 


1.310 


± 


0.090 


1.240 


± 


0.090 


6400 ± 


100 


WASP-4 


1 


,33823 


0.02255 


1.602 


+ 


0.258 


1.237 ± 


0.064 


0.873 


+ 


0.036 


0.850 


+ 


0.120 


5500 ± 


100 


WASP-5 


1 


.62843 


0.02670 


1.100 


± 


0.070 


1.580 ± 


0.110 


1.030 


+ 


0.060 


0.970 


+ 


0.090 


5700 + 


150 


WASP-6 


3 


,36101 


0.04210 


1.220 


± 


0.050 


0.500 ± 


0.030 


0.870 


+ 


0.030 


0.880 


± 


0.070 


5450 ± 


100 


WASP-8 


8 


,15875 


0.07930 


1.170 


+ 


0.100 


2.230 ± 


0.100 














5600 ± 


300 


WASP- 10 


3 


,09276 


0.03781 


1.080 


± 


0.020 


3.150 ± 


0.120 


0.698 


± 


0.012 


0.750 


± 


0.030 


4675 ± 


100 


WASP- 12 


1 


,09142 


0.02290 


1.790 


+ 


0.090 


1.410 ± 


0.090 


1.570 


+ 


0.070 


1.350 


+ 


0.140 


6250 ± 


150 


WASP- 13 


4 


,35298 


0.05270 


1.220 


± 


0.130 


0.460 ± 


0.060 


1.340 


± 


0.130 


1.030 


± 


0.110 


5826 ± 


100 


WASP- 15 


3 


.75207 


0.04990 


1.430 


± 


0.080 


0.540 ± 


0.050 


1.477 


± 


0.070 


1.180 


± 


0.120 


6300 ± 


100 


WASP- 16 


3 


,11806 


0.04210 


1.019 


+ 


0.072 


0.855 ± 


0.045 


0.946 


+ 


0.060 


1.020 


+ 


0.070 


5700 ± 


150 


WASP- 18 





.94145 


0.02026 


1.115 


± 


0.063 


10.300 i 


: 0.690 


1.216 


± 


0.070 


1.250 


± 


0.130 


6400 ± 


100 


CoRoT-2 


1 


.74300 


0.02800 


1.465 


+ 


0.029 


3.310 ± 


0.160 


0.902 


+ 


0.018 


0.970 


+ 


0.060 


5625 ± 


120 


CoRoT-3 


4 


.25680 


0.05694 


1.060 


± 


0.120 


21.660J 


: 1.000 


1.540 


± 


0.090 


1.360 


± 


0.090 


6740 ± 


140 


CoRoT-4 


9 


.20205 


0.09000 


1.190 


± 


0.060 


0.720 ± 


0.080 


1.170 


± 


0.020 


1.160 


± 


0.030 


6190 + 


: 60 


CoRoT-5 


4 


,03790 


0.04947 


1.388 


± 


0.046 


0.470 ± 


0.050 


1.186 


± 


0.040 


1.000 


+ 


0.020 


6100 + 


: 65 


OGLE-TR-L9 


2 


,48553 


0.03080 


1.610 


± 


0.040 


4.500 ± 


1.500 


1.530 


+ 


0.040 


1.520 


± 


0.080 


6933 + 


:60 


OGLE-TR-10 


3 


,10128 


0.04162 


1.245 


+ 


0.095 


0.610 ± 


0.130 


1.140 


+ 


0.080 


1.100 


+ 


0.050 


6075 + 


:86 


OGLE-TR-56 


1 


,21191 


0.02250 


1.300 


± 


0.050 


1.290 ± 


0.120 


1.320 


± 


0.060 


1.170 


± 


0.040 


6119 + 


: 62 


OGLE-TR-111 


4 


,01445 


0.04670 


1.010 


+ 


0.040 


0.520 ± 


0.130 


0.831 


+ 


0.031 


0.810 


+ 


0.020 


5044 + 


: 83 


OGLE-TR-113 


1 


,43248 


0.02290 


1.090 


+ 


0.030 


1.350 + 


0.190 


0.770 


± 


0.020 


0.780 


+ 


0.020 


4804 ± 


106 


OGLE-TR-132 


1 


,68987 


0.02990 


1.180 


+ 


0.070 


1.140 ± 


0.120 


1.340 


± 


0.080 


1.260 


± 


0.030 


6210 + 


: 59 


OGLE-TR-182 


3 


,97910 


0.05100 


1.210 


± 


0.160 


1.010 ± 


0.115 


1.140 


± 


0.150 


1.140 


± 


0.050 


5924 + 


: 64 


OGLE-TR-211 


3 


,67724 


0.05100 


1.395 


± 


0.125 


1.030 ± 


0.200 


1.640 


± 


0.130 


1.330 


± 


0.050 


6325 + 


: 91 



Table 2. Parameters of the considered transiting planetary systems. 



where bo and cq are free coefficients, J -3/2 and J3/2 are Bessel 
functions of the first kind of order -3/2 and 3/2, respectively, 
q = \a\r, and qo = \a\R. Making use of Eq. ©, the magnetic 
field components are: 

R 2 

B r = 2B -rg(q)cos9, 



R 2 

B = -B Q \a\—g'(q)sm6, (5) 
r 

R 2 

B^ = aB — g{q) sin 0, 
r 

where g'(q) = dg/dq. A linear force-free field as given by 
Eqs. (0 extends to the infinity with an infinity energy. We con- 
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Name 


v sin i 


^rot Age 


^synch 


e 


References 




(kms 


1 . 

') 


(days) (Gyr) 


(Gyr) 






HD 149026 


6.00 ± 


0.50 




2.00- 3.60 


34.10 


0.0 


Wo07 Sa05 


HD 189733 






11.953 + 


: 0.009 




0.004 ± 0.003 


He07 Mo07 Po07 Wi07 


HD197286/WASP7 


17.00 ± 


: 2.00 








0.0 


He09 


HD209458 






11.4 + 


: 1.5 




0.015 ± 0.005 


Wi05 Wi05a 


TrES-1 


1.30 ± 


0.30 








0.0 


Na07 Wi07a A104 


TrES-2 


1.00 ± 


0.60 




2.80- 7.80 


0.56 


0.0 


Wi08 So07 


TrES-3 


1.50 ± 


1.00 




0.10- 3.70 


0.04 


0.0 


So09 


TrES-4 


8.50 ± 


0.50 




2.50- 4.40 


10.45 


0.0 


So09 


XO-1 


1.11 ± 


0.70 




0.60- 5.50 


9.31 


0.0 


Ho06 Mc06 


XO-2 


1.30 ± 


0.50 








0.0 


Fe09 Bu07 


XO-4 


8.80 ± 


0.50 




1.50- 2.70 


12.03 


0.0 


Mc08 


XO-5 


0.70 ± 


0.50 




12.80-16.80 


1.88 


0.010 ± 0.013 


Pa09 


HAT-P-1 


3.75 ± 


0.58 




1.60- 4.60 


130.20 


0.0 


Jo08 Wi07a Ba07a 


HAT-P-3 


0.50 ± 


0.50 




0.10- 6.90 




0.0 


To07 


HAT-P-4 


5.50 ± 


0.50 




3.60- 6.80 


7.87 


0.0 


To08 Ko07 


HAT-P-5 


2.60 ± 


1.50 




0.80- 4.40 


2.53 


0.0 


To08 Ba07b 


HAT-P-6 


8.70 ± 


1.00 




1.60- 2.80 


28.17 


0.0 


To08 No08 


HAT-P-8 


11.50 ± 


: 0.50 




2.40- 4.40 


7.33 


0.0 


La08 


HAT-P-9 


11.90 ± 


: 1.00 




0.20- 3.40 


204.30 


0.0 


Am09 Sh09 


HAT-P-10/WASP11 


0.50 ± 


0.20 




4.10-11.70 


11.05 


0.0 


We09 Ba09a 


HAT-P-1 2 


0.50 ± 


0.40 




0.50- 4.50 


38.73 


0.0 


Ha09 


HAT-P-1 3 


4.10 ± 


0.50 




4.20- 7.50 


2.79 


0.021 ± 0.009 


Ba09b 


WASP-1 


5.79 ± 


0.35 




1.00-3.00 


3.11 


0.0 


Ca07 St07 Sh07 Ch07 


WASP-2 












0.0 


Ca07 St07 Sh07 Ch07 


WASP-3 


13.40 ± 


: 1.50 








0.0 


Gi08 Po08 


WASP-4 


2.00 ± 


1.00 








0.0 


Gi09 


WASP-5 


3.50 ± 


1.00 








0.038 ± 0.026 


Gi09 


WASP-6 


1.60 ± 


0.30 








0.054 + 0.018 


Gi09a 


WASP-8 












0.0 


Sm09 


WASP- 10 






11.91 + 


: 0.05 




0.059 ± 0.014 


Jo09 Sm09 


WASP- 12 


2.20 ± 


1.50 








0.049 ±0.015 


He09 


WASP- 13 


2.50 ± 


2.50 




3.60-14.00 




0.0 


Sk09 


WASP- 15 


4.00 ± 


2.00 




2.60- 6.70 


22.24 


0.0 


We09a 


WASP- 16 


3.00 ± 


1.00 




0.10- 8.10 


12.21 


0.0 


Li09 


WASP- 18 


11.00 ± 


: 1.50 




0.50- 1.50 


6.5 x 10~ 4 


0.009 ± 0.003 


He09a 


CoRoT-2 






4.52 ± 


0.15 0.10- 1.70 


0.36 


0.0 


Bo08 A108 


CoRoT-3 


17.00 ± 


: 1.00 




1.60- 2.80 


1.41 


0.008 ±0.015 


De08 Tr09 


CoRoT-4 






9.20 + 


: 0.3 0.70- 2.00 


2728.00 


0.0 ± 0.100 


Ai08 Mo08 La09 


CoRoT-5 


1.00 + 


1.00 




5.50- 8.50 




0.090 ± 0.090 


Ra09 


OGLE-TR-L9 


39.33 i 


:0.40 




0.10- 0.70 


0.92 


0.0 


Sn09 


OGLE-TR-10 


7.70 ± 


3.00 








0.0 


To04 Bo05 Ko05 Po07a 


OGLE-TR-56 


3.00 ± 


3.00 








0.0 


To04 Bo05 Ko05 Po07a 


OGLE-TR-111 












0.0 


Po04 Mi07 


OGLE-TR-113 


3.00 ± 


3.00 








0.0 


Bo04 Di07 Ko04 


OGLE-TR-132 


3.00 ± 


3.00 




0.10- 1.40 




0.0 


Gi07 Mo04 Bo04 


OGLE-TR-182 












0.0 


Po08a 


OGLE-TR-211 












0.0 


Ud08 



Table 3. Parameters of the considered transiting planetary systems. 



Reference codes: Ai08: lAigrain et all (2003) : A104: lAlonso et alj t25f3l : MPS: IMonso et all |200H) : Am09: [Ammler-von Eiffel al. 2009): Ba07: Bakos et al. 20073); Ba07a: IBakos et al. 2007bj); Ba07b: 
IBakos el all (20073) : Ba09a: lBakoFeTall 12009a ) : Ba09b: Bakos et ail l20O9U) : Bo04: IBouchveTai1 12001 ; Bo05: IBouchv"eTai] 12003) : Bo08: IBouchfetai]l200l : Bu07: IBurke eTaTl (2TO7I) : CaQ7: ICameron et all 

(2001 : Ch07: ICharbonneau el all (20071 ); De08: IDeleuiieTai] (2001 ; Di07: IDi'az et"ai] (20071 ); Fe09: 1 Fernandez el all [20091) ; Gi07: IGillon etaTl (20Ol ; Gi08: IGibson'eta!l 12001 ; Gi09: IGiHoneTail 12009 J ; 

Gi09a:lGUlon et al. 2009b}); Ha09: IHartman et all 12009 J) ; He08: IHenrv & Winnl 12001 ; He09: IHebb eTaTl 120091) ; He09a: lHeifier el all I2009d) ; He09b: lHemer et all I2009U) : Ho06: IHolman et all 120061) ; Jo08: 
IJohnsonet all (20081) : Jo09: IJoshi et ail (20091) ; Ko04:lKonacki et al. 1 2004): Ko05: Konad JeTaTl 155531 ; Ko07: IKovacs"etail (20071) ; La08: ILalham"etail 120081) ; La09: ILanza etai] (SfjOll : Li09 : l Lister et"all 120091) ; 

Mc06: IMcCullouehetall (2001 ; Mc08:[McCullough et 3] (20081 ); Mi07: IMinniti7Tai] l20O7t ; Mo04:[Moutou et al. 1 200§; Mo07: IMoutou et all (2001 ; Mo08: IMoutou et all (20081 ; Na07: INarita eTail (2001 ; 

No08: INoves"eTall 120081) ; Pa09:|Pal et al. 2009): Po04: Pont et ail |2003| : Po07: IPont etaT] I2007al) : Po07a: |PonTeTail I2007H ; Po08: IPollacco et all (20081) : Po08a: |Ponn=i"aTl (20081) ; Ra09: IRauer eTail 12009]) ; 

Sa05: ISatoetan(200l ; Sh07: IShporererail (2001 : ShQ9: IShDoreretail 12001 ; Sk09: ISkillen eTail (20091) : Sm09: ISmith eTail 12009]) ; Sn09: ISnellerreTai]l2SSl : So07: ISozzetti et all (2001) : So09: ISozzettT"eTan 

[2009] ); St07: IStempels et all (5)071 ); To04: ITorres eTail (2001 ; ToQ7: ITorres eTaTl 120071) ; To08: ITorres eTail (20081 ); TrQ9: ITriaud etail j200Sh ; Ud08: IUdalski et all (20081 ); We09: IWest eTaTl I2009J) ; We09a: 
IWest et all I2009M) : Wi05: IWinn eTail 12003) ; Wi05a: rwittenmver et all (2001 ; Wi07: IWinn eTail I2007tl) ; Wi07a: IWu7neTail I2007M : Wi07b:|Wmet2] I2007d) ; Wi08: IWinn eTail I2008al) ; Wo07: 1 Wolf eTail 
2007). 
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sider its restriction to the radial domain < q < q^, where q^ 
is the first zero of g{q), s o that all the magneti c field lines are 
closed in our model (see Chandrasek harl 1 9561 for the bound- 
ary conditions at r — = q^/\a\). 

The magnetic field geometry specified by Eqs. ((5) depends 
on two independent parameters, i.e., a and bo/co. They can 
be derived from the boundary conditions at the stellar photo- 
sphere, i.e., knowing the magnetic field B (s, (Q, <f>) on the sur- 
face at r — R. Using the orthog onality properties of the basic 
poloidal and toroidal fields (see IChandrasek har 1961), we find: 



-rB () R 2 



— \a\B R 3 g'(q ) 



-aBoR 



= f *? 

--£ 

= - r < 



cos QdL, 



sinQdL, 



sin QdL, 



(6) 



where L(R) is the spherical surface of radius R, and dL = 
R 2 sin QdQd(f>. Note that the photospheric magnetic field compo- 
nents can be measured by means of spectropolarimetric tech- 
niques if the star rotates fast enough (v sin i > 10 — 15 km s _1 ) 
as shown in t he cas e of, e.g., t Boo bv lCatala et al. (2007) and 
iDonati et al.l d2008l) . With such a kind of observations, Eqs. © 
can be applied to derive the parameters of the coronal field 
model and its topology, provided that the field is approximately 
force-free down to the stellar photosphere. We shall refer to 
this approach in Sect. 14.21 where our model will be used to 
study the evolution of stellar angular momentum. 

The magnetic energy E of the field confined between the 
spherical surfaces r — R and r — can be found from Eq. (79) 
in § 40 of IChandrasekharl d 1 96 ll) : 



E = E p [2 + qmiXgXqOf ~ qfaXqo)? ~ <?o} , (7) 

where E p = (47r/3fi)R 3 B 2 l is the energy of the potential dipole 
field with the same radial component at the surface r — R, and 
/i is the magnetic permeability . The relative magnetic helicity 
//r, as defined bv lBergei ( 1985 ). can be found from his Eq. (19) 
and is: 



V(9o) + 



8tt E 



3 qoEp 



W 
a 



(8) 



Note that the field obtained by changing the sign of a has the 
same poloidal components B r and Bg, and energy E, while the 
toroidal component B^ and the relative helicity Hr become op- 
posite. Further inform ation on the field described by Eqs. (0 
can be found in, e.g., lLanzal (12009b . 

For a finite a, E > E p because the potential field has the 
minimum energy for a given Br . If we consider all magnetic 
fields with one end of their field lines anchored at r = R and 
the other out to the infinity, satisfying the same boundary con- 
ditions of our field at r — R, the field with the lowest possible 
energy is called t he Aly field and its energy E\\ y = 1 .66£p (see 



Fiver et al. 2004). We assume that the Aly energy is an upper 
bound for the energy of our field because it is the lowest energy 
allowing the field to open up all its lines of force out to the in- 
finity driving a plasma outflow similar to a solar coronal mass 
ejection. 




Fig. 6. Upper panel: The outer radius of the coronal field 
at the Aly energy limit vs. the absolute value of the force-free 
parameter a; middle panel: the absolute value of the relative 
magnetic helicity of the field at the Aly energy limit \Hr\, as 
computed from Eq. (|S), vs. \a\; lower panel: the value of the 
parameter bo/co corresponding to the Aly energy limit vs. \a\. 



Since all magnetic field lines are closed in our model, it is 
not possible to have a steady flux of angular momentum toward 
the infinity as, e.g., in the open field configuration of solar coro- 
nal holes. Therefore, we assume that the loss of angular mo- 
mentum occurs only when the field energy reaches the Aly en- 
ergy and the coronal field opens up toward the infinity driving 
out all the coronal mass. The rate of angular momentum loss 
depends on the angular momentum stored in the coronal field 
and the rate of occurrence of such events that we call coronal 
mass ejections (CMEs) by analogy with similar solar events 
which usually involve only a single active region and not the 
whole coronal field as in our simplified model. 

If we fix the value of the parameter a, the value of the ra- 
tio bo/ Co corresponding to the Aly energy can be determined 
numerically. We plot in Fig. [6] the outer radius, the relative 
magnetic helicity and the value of bo/co vs. a for a field at the 
Aly energy limit. There is a remarkable decrease of the outer 
field radius and of the relative magnetic helicity with increas- 
ing a. This implies that a coronal configuration with a greater 
a is more tightly confined than one with a lower value of a. 
Moreover, its relative helicity is significantly smaller than in 
the case with a lower a. 

3.2. Some considerations on non-linear force-free field 
models 

The application of our linear model is justified in view of its 
mathematical simplicity and the hypothesis of a corona with 
closed field lines. A general treatment of the coronal field topol- 
ogy poses formidable mathematical problems even consider- 
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ing only force-free fields. Therefore, we restrict ourselves to 
a special class of non-linear force-free models to make some 
progress. 

Low & Lou^OOOh . lFiver et all J^OOl . lZhang et al](l2006l) . 
and other authors, showed that a force-free field can extend 
to the infinity with a finite magnetic energy and relative he- 
licity if a is not uniform, i.e., it is a non-linear force-free 
field. Such fields are suitable to model a stellar corona without 
the limitations of our adopted linear force-free model, specif- 
ically the property that all field lines are closed. An example 
of such a non-line ar axisymmetric field has been provided by 
Low & Lou (2000). It is obtained by assuming Q(A) = OA « 
in Eq. (0, with being a constant and A(r, 9) = P(ff)r~"; i.e., 
the flux function A is of separable form, with P being a con- 
tinuous function of 9 defined in the closed interval [0, n] which 
must vanish at the poles to ensure the regularity of the field. 

Let us consider a magnetic field line going from the sur- 
face of the star to the infinity. Since A must be constant on it, 
we conclude that P must be zero, otherwise P would increase 
without bound following the field line for n > 0, which is of 
course not possible in view of its continuity. In other words, all 
the magnetic field lines going out to the infinity must be rooted 
on the surface of the star at colatitudes where P(9) = 0. 

Isolated field lines going out to the infinity do not have 
any effect on the stellar angular momentum loss because the 
mass flow along them vanishes. To produce an angular mo- 
mentum loss, we need a flux tube with a finite cross-section, 
i.e., a bundle of field lines covering a finite colatitude interval 
on the surface of the star, say, 9\ < 9 < 92, in which P(9) = 0. 
Considering the expression for B r in Eq. (0, it follows that 
B r = inside that interval because (dA/d9) r= R = 0. In other 
words, those field lines do not intersect the surface of the star 
which is the source of any wind mass loss. We deduce that there 
can be no steady mass loss along those field lines, hence their 
contribution to the angular momentum loss is negligible. 

In conclusion, a non-linear force-free field with a separa- 
ble flux function of the kind proposed by iLow & Lou (2000), 
although endowed with field lines going from the surface of 
the star to the infinity, cannot sustain a steady angular momen- 
tum loss from the star through a continuous wind flow. The 
only available mechanism is therefore represented by the CME 
events considered above. 

Of course such a conclusion has been obtained for a spe- 
cific class of non-linear models and cannot be generalized to 
all possible non-linear fields. Nevertheless, it shows that there 
are non-linear force-free fields for which our approach of con- 
sidering only the angular momentum stored into a closed-field 
corona is valid. Our linear model has the advantage of a spher- 
ically symmetric outer boundary of the corona at r — r^, while 
non-linear models may have a boundary that depends on the 
colatitude. 



3.3. Angular momentum content of the coronal field 

To compute the angular momentum content of the plasma 
trapped into the closed coronal field introduced in Sect. 13.11 
we need to evaluate its moment of inertia. For the sake of sim- 



plicity, let us assume that the plasma has a uniform temperature 
T and is in hydrostatic equilibrium along each magnetic field 
line. The eq uation of hydrostatic equilibrium reads (cf., e.g., 
|Priestlll982l) : 



s ■ (-Vp + pVO) = 0, 



(9) 



where s is the unit vector in the direction of the magnetic 
field, p the plasma pressure, p its density and €> the total po- 
tential including the gravitational and centrifugal terms. Since 
p = (R/p)pT and <D = GM/r + 5Qr 2 sin 2 6», where R is the 
gas constant, p the molecular weight of the plasma, G the grav- 
itation constant, and M the mass of the star, Eq. (0 can be 
immediately integrated to give: 



p = po exp 



^(D-cDo) 



(10) 



where po and <t>o are the density and the potential on a field line 
at the base of the corona, i.e., at a radius r = ro. The density 
is maximum on the equatorial plane of the star, i.e., at 9 — n/2 
because there €> is maximum. For our purposes, we can esti- 
mate the maximum moment of inertia, which gives us an upper 
limit for the angular momentum loss, by assuming a spherically 
symmetric distribution of the density computed by evaluating 
O on the equatorial plane, i.e.: 



/>!/■) =/>oC\p^ 



(?-) 



+ £ 



( r 2 
V0 



(11) 



where H p = RT/(figo) is the pressure scale height at the base of 
the corona, with g = GMj r 2 , being the gravitational accelera- 
tion at r — ro, and e = r^Q. 2 /(2GM) the ratio of the centrifugal 
to the gravitational potential at the base of the corona. We as- 
sume that the base of the corona is at ro = 2R, outward of which 
the plasma temperature is assumed to be constant. Considering 
a star with the radius and mass of the Sun and a mean molecu- 
lar weight p = 0.6, we have H p a 2.04 x 10 8 (r/10 6 ) m, with T 
in K. As a typical coronal density, we adopt the mean ele ctron 



10 1Z nr 3 dCoxl 



density of the solar corona at 2R Q , i.e., n 
20001) . 

The moment of inertia of the coronal plasma can be easily 
computed noting that in our assumptions the density is a func- 
tion of r only. If the corona extends up to the limit radius r^, its 
moment of inertia is: 



= \ pr 2 sin 2 9 dV, 
Jv 



(12) 



where V is the volume of the corona, that is the spherical shell 
between radii ro and r^. 

In Fig. [7] we plot the moment of inertia of a corona at 
the Aly energy limit vs. the force-free parameter a for a star 
analogous to the Sun, setting the coronal base at ro = 2R with 
« e = 10 12 rrT 3 . We compute / for two temperature values, i.e., 
T = 1.6 x 10 6 K, typic al of stars with a low level of coronal 
emission as the Sun (cf. ICoxll2000h . and T = 3 x 10 6 K, which 
is characteristic of stars with a moderately high level of coronal 
emission, i.e., with an X-ray flux about one order of magnitude 
great er than the Su n at the maximum of the 11 -yr cycle (cf., 
e.g., ISchmittlll997h . We assume that the rotation period P r ot is 
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Fig. 7. The moment of inertia of the coronal plasma vs. the ab- 
solute value of the force-free parameter a, computed accord- 
ing to Eqs. (Hi} and (|T2j for a star with R - R & , M - M B , 
ro = 2R Q , a base electron density n e (ro) = 10 12 rrT 3 , and tem- 
perature T = 1.6 x 10 6 K (solid line) or T = 3 x 10 6 K (dotted 
line). The linear force-free configuration of the coronal field is 
that corresponding to the Aly energy limit for the given value 
of a. The slope changes at |a| = 0.2 are due to the fixed value 
ofP rot = 3 days for M > 0.2. 



inversely correlated with |or|, as will be discussed in Sect. 14.21 
Specifically, we assume that P rot increases linearly between 3 
and 24 days when \a\ decreases from 0.2 to 0.025. 

The plots in Fig. [7] are terminated where the potential en- 
ergy or the internal energy of the plasma exceed 0. 1 of the total 
magnetic energy of the field computed for Z?o = 20 G because 
the force-free condition is no longer valid in such a case. The 
potential energy Ecc and the internal energy U of the coronal 
plasma are evaluated as: 



Ecc 



and 



U 



f p<SdV, 
Jv 



1 RT 



7- 1 P- Jv 



pdV, 



(13) 



(14) 



where y — 5/3 is the ratio of the specific heats of the plasma. 

Note the decrease of the moment of inertia by ~ 2.5 or- 
ders of magnitude when \a\ increases from 0.06 to 0.2 owing to 
a remarkable decrease of the outer radius of the corona (cf. 
Fig. [6] upper panel). On the other hand, an increase of the coro- 
nal temperature by a factor of ~ 2, as expected for rapidly ro- 
tating stars, produces an increase of the moment of inertia only 
by a factor of ~ 3 — 4. The moment of inertia is directly pro- 
portial to the base density, so a change of po by, say, one order 
of magnitude produces a corresponding change in the moment 
of inertia. We conclude that the most relevant variation of the 
moment of inertia of the corona is produced by a variation of 
the force-free parameter a. 



4. Application to stellar angular momentum 
evolution 

4.1. Pre-main-sequence evolution 



To apply the results of Sect. I3.3I to the problem of stellar an- 
gular momentum evolution, we need to define the initial rota- 
tion state of a star. Stars with hot Jupiters are accompanied by 
circumstellar discs during the pre-main-sequence (hereinafter 
PMS) phase of their evolution, which play a fundamental role 
in the formation and orbital evolution of their planets. The an- 
gular velocity of a PMS star is equal to the Keplerian angular 
velocity of its disc at the so-called corotation radius. It is lo- 
cated ~ 5 - 10 percent outside the inner boundary of the disc, 
where it is truncated by the stellar magnetic field (cf., e.g.. 



Cameron & Campbell 1 119931: iTinker et a l. 2002; Scholz et al 



2007; iBouvierf 20081 and references therein). According to 



the current theoretical scenario, hot Jupiters are formed at 
several AUs from their stars, beyond the snow line where 
volatile elements can condense, and then migrate toward their 
stars on a timescale not exc eeding 10 5 - 10 6 yr (cf., e.g., 
Papaloizou & Terquem 120061) . If the stellar magnetic field is 



strong enough, the Keplerian shear induces a sizeable toroidal 
magnetic field in the ionized region of the disc close to the 
star and the inward migration of the plane t may be halte d 
close to the corotation radius, as suggested bv lTerqueml d2003h . 
The field intensity required to halt inward migration depends 
on the variation of the parameter ft vs. the radius within the 
disc, where ft = (c s / va) 2 with c s bein g the sound speed and 
va the Alfven speed. Terqueml (2003) showed that values of 
/3 ~ 10-100 can be sufficient to halt planetary migration. 

According to this scenario, the initial rotation period of 
the star P rot is approximately equal to the orbital period of 
the planet, i.e., it is between 3 and 10 days. The typical life- 
time of the disc does not exceeds 5-10 Myr which is shorter 
than the timescale of contraction to reach the zero-age main se- 
quence ( hereinafter ZAMS) for stars having a mass lower tha n 
* 2M^ dTinker et alJ l2002t iBouvierl \20oA lMamaiekll2009h . 
When the disc disappers, stellar rotation is no more locked 
and the rotation period decreases during the approach to the 
ZA MS owing to the reduction of the moment o f inertia of the 
star dScholz et alj|2007t llrwin & Bouvierj|2009l) . In Fig. |8] we 
plot the evolution of the radius a nd the moment of inertia dur- 
ing the PMS phase, according to lSiess et al . (2000), for stars of 
1 .0, 1 .2, and 1 .4 M , respectively. The decrease of the moment 
of inertia is greater than expected on the basis of the contraction 
of the radius because the internal structure changes also with an 
increase of the mass of the radiative core as the star approaches 
the ZAMS. The reduction of the moment of inertia occurring 
between disc decoupling and arrival onto the ZAMS is by a fac- 
tor of ~ 5 if the disc lifetime is 5 Myr. This implies a remark- 
able acceleration of stellar rotation which destroys any syn- 
chronization with the planetary orbit attained during the pre- 
vious disc-locking phase. This would give «/Q ~ 0.2 - 0.3 for 
ZAMS sun-like stars with hot Jupiters, for which case there is 
no evidence in our sample of transiting planets. Therefore, we 
conjecture that some process is at work to restore synchroniza- 
tion when a planet-harbouring star is approaching the ZAMS. 
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A candidate m echanism is a magne tocentrifugal stellar wind, as 
suggested bv Lovelace et al.1 (120081) . Considering a star which 
was released by its disc with a rotation period of 8 days, it 
would reach the ZAMS with a period of only 1.6 days, if the 
reduction of the moment of inertia is not counteracted by any 
other process. Assuming that the young contracting star has a 
surface magnetic field of 10 3 G, the torque exerted by its coro- 
nal field on the planet would transfer most of the stellar angular 
momentum to the planet itself on a time scale of 3 - 5 Myr. 



restori ng a synchronous rotation state. Recently. IVidotto et al 



(2009) have revisited such a mechanism considering a more 
re alistic wind model tha n the Weber & Davis model adopted 
bv lLovelace et al. (2008). They find timescales longer by one 
order of magnitude for the angular momentum exchange be- 
tween the star and the planet, which are still acceptable in the 
framework of our model. We conclude that a magnetocentrifu- 
gal wind may maintain synchronization in solar-like stars ac- 
companied by a hot Jupiter during PMS evolution after the star 
has been released by its disc. This implies that the star arrives 
on the ZAMS in an approximate synchronous state of rota- 
tion. After the star has settled on the ZAMS, the efficiency of 
the stellar hydromagnetic dynamo decreases with respect to its 
PMS phase because the volume of the outer convection zone 
is significantly smaller than in the PMS phase, so the magnetic 
field intensity at the surface drops and the coupling provided 
by the magnetocentrifugal wind virtually vanishes. From this 
point on, the evolution of the spin and the orbital angular mo- 
mentum are decoupled and we can study the evolution of stellar 
rotation treating the angular momentum loss from the corona 
by means of the model of Sect. 13.31 



In addition to the scenario proposed above, another evolu- 
tionary sequence is possible if the magnetic field of the star 
truncates the disc and couples the rotation of the star to its in- 
ner edge, but it is not strong enough to halt the migration of the 
planet (i.e., f3 > 10 - 100). In this case, the planet will continue 
to migrate inward until its orbital period becomes half of the 
period at the corotation radius because the angular momentum 
exchange betwe en the planet a nd the disc proceeds via the 2:1 
resonance (see iLin et al. I l 19961 and references therein). In this 
case, the initial rotation period of the star is twice the orbital 
period of its hot Jupiter. If the star is massive enough, say at 
least M x 1.5 M , and its disc is long-lived, say, * 10-15 Myr 
(Mamaie ket alJl2002l) . it can reach the ZAMS while still be- 
ing locked to its disc, thus starting its evolution in a rotational 
status with n/Q 2. If those stars do not appreciably loose 
angular momentum during their main-sequence evolution (see 
Sect. 14.2b . this may explain the observed concentration of sys- 
tems with r e ff ^ 6200 K around n/Q = 2 seen in Figs.[T]and[2] 



4.2. Main-sequence evolution 

Starting from the initial status on the ZAMS described in 
Sect. 14.11 we want to account for the main features of the n/Q 
distribution found in Sect. [2] namely the dependence of n/Q 
on the effective temperature, with stars having T e ff ^ 6000 K 
showing a generally smaller n/Q than cooler stars, and the de- 




10 20 30 

Age (Myr) 

Fig. 8. Upper panel: Radius of PMS stars of different mass vs. 
time measured from their birth line; the radius is normalized 
to the value at an age of 5 Myr, corresponding to the average 
lifetime of the circumstellar discs; different linestyles refer to 
different masses: solid: 1 M ; dotted: 1.2 M G ; dashed: 1.4 M . 
Lower panel: Moment of inertia of PMS stars vs. time from 
their birth; the moment of inertia is normalized at the value at 
an age of 5 Myr; different linestyles refer to different masses as 
in the upper panel. 



pendence of n/Q on the stellar rotation period, found in stars 
having r eff > 6000 K. 

To study the angular momentum evolution on the main 
sequence, we apply the model of Sect. 13.31 The force-free 
paramete r a of the co ronal field plays a crucial role in that 
model. lLanzal (120081) proposed a method to estimate a in 



stars showing chromosp heric hot spots rotating sync hronously 
with their hot Jupiters dShkolnik et al] 120051 Hooi) . To date, 
only five stars have been modelled, so conclusions based on 
such a method are still preliminary. Nevertheless, for F-type 
stars, i.e., HD 179949, v Andromeda e, and t Boo, havin g 
r eff > 6200 K and P mt < 12 days (cf. IShkolnik et a"i1l2008l) . 
the values of \a\ fall between 0. 1 and 0.2, while for the two K- 
type stars HD 189733 and HD 192263, ha ving T eS * 5000 K 



and P rot > 12 days (cf. ISantos et al ] l2003l) . \a\ ranges between 



0.025 and 0.1. Note tha t such values we re obtained with the 
non-force-free model of Neukirch] ( 1995 ). but the typical val- 
ues of a obtained with a purely force-free model do not differ 
by more that 10-20 percent. 

A motivation for a greater value of \a\ in F stars than in G 
and K stars may be the stronger toroidal field at their surface 
produced by a greater relative differential rotation. Assuming 
that the measured photospheric field is a good proxy for the 
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field at the base of the corona, an estimate of a can be obtained 
by comparing the first and the third of Eqs. (O which yields: 



a = — 



(15) 



Indeed spe ct ropola rimet ric observa t ions of t Boo by 
Donati et all d2008b and iFares et al.l d2009b show that the 



star has an oscillating field with a predominantly toroidal com- 
ponent during a significant fraction of its activity cycle. This 
may be a good example of a fast rotating F-type star with a hot 
Jupiter because its mean rotation period is 3.3 days and it also 
shows a surface differential rotation with a relative amplitud e 
of ~ 0.2 between the equator and the pole (Cat ala et alJ l2007). 
In the case of CoRoT-4a, time series spot modelling suggests a 
surface differential rotation comparable to that of t Boo, thus 
supp orting the presenc e of a predominantly toroidal surface 
field dLanza et al.ll200S " 



Barn es et ai](l2005h and iReinersI J2006I) show that the 
amplitude AQ of the surface differential rotation decreases 
strongly with the decrease of the effective temperature of the 
star, viz. AO oc r 8 „ 9±0 - 3 . Recent spectropolarimetric observa- 



tions bv lPetit et al.l (120081) indicate that the photospheric mag- 
netic field of G-type stars (r eff - 5700 - 6000 K) with a ro- 
tation period below 10-12 days is predominantly toroidal, 
while stars with a longer rotation period have a predominantly 
poloidal field. Such a contrast may come from a different am- 
plitude of the shear at the boundary between the radiative 
core and the convective envelope (also called the tachocline 
in the Sun), which m ay be greater in hotter and fast-rotating 



stars. iBouvi er (2008), specifically considering stars with mas- 
sive planets, suggested that their lower Lithium abundance may 
be the result of an enhancement of the turbulence at the core- 
envelope interface induced by hydrodynamic or magnetohy- 
drodynamic instabilities associated to a sizeable shear local- 
ized at the interface. This suggests that a sizeable toroidal field 
is present in such stars, at least during the first phase of their 
evolution on the main sequence, produced by the shearing of a 
radial poloidal field close to the base of their convection zones. 

Considering rapidly rotating (P lot < 10-12 days) F-type 
stars, we adopt 0.1 < \a\ < 0.2, yielding a typical moment of 



inertia of their coronae ranging from ~ 6 x 10 to ~ 10 kg m 
(cf. Fig. |7). The timescale for angular momentum loss can be 
estimated as: 



/„ y-MR 1 
taml = yAf = — J — At > 



(16) 



where /» = y 2 MR 2 is the moment of inertia of the star, yR ~ 
0.35/? is its gyration radius, and Af is the mean time interval 
between the CME events that produce the loss of the angular 
momentum of the stellar corona. 

We can estimate a lower limit for Af from the ratio between 
the total energy of the coronal field and the X-ray luminosity of 
the star, i.e., At a E/L x . For r Boo and HD 179949, the average 
value of L x is ~ 3 x 10 21 W, while E = 7.6 x 10 27 J at the Aly 
limit for Bo = 10 G (a field intens ity measured in r Boo by 



Donati et al]|2008t IFares et al.ll2009ft . giving Af ~ 2.6 x 10 6 s. 



global coronal field topology, as discussed in lLanza (2009), i.e., 
Af ~ 2.6 x 10 7 s, or ~ 300 days. We shall adopt Af = 2.6 X 10 6 s 
because a variation of Af in Eq. dT6l l can be compensated by a 
change of / given that the density po at the base of the corona 
may vary by one order of magnitude. 

In the case of r Boo, R = 1.6 Rg, M = 1.2 M , and |or| = 
0.12, yielding / ~ 3 x 10 35 kg m 2 ; thus we find an angular mo- 
mentum loss timescale of the order of taml ~ 100 Gyr. Such 
a value implies that the initial angular momentum of r Boo 
remains approximately constant during its main-sequence life- 
time. In other words, the observed synchronization between 
the average stellar rotation and the orbital period of the planet 
should be a remnant of the initial state of the system when the 
star settle d on the ZAMS. A similar conclusion is reached for 
CoRoT-4 dLanza et alj20 09). An angular momentum loss time 
scale of the order of 100 Gyr accounts also for the rotation peri- 
ods of the mid-F type stars in the systems XO-4 and HAT-P-6, 
which again appear to be remnants of their ZAMS rotational 
status, in this case with an initial n/Q ^ 2. 



In the light of the results of lPetit etall d2008l) . stars with 



' eft' 



An upper limit may come from the timescale for changing the 



6000 K and rotation periods longer than ~ 10 days 
should be characterized by a smaller surface toroidal field than 
more rapidly rotating stars which implies a smaller value of a. 
Therefore, their angular momentum loss time scale is expected 
to be shorter than that of the rapidly rotating F-type stars con- 
sidered above, which may account for the dispersion of n/Q 
observed in the effective temperature range 6000 - 6500 K. For 
systems such as WASP-1, WASP- 15, or WASP- 18, adopting 
R = 1.4 R©, M = 1.2 Mq, and \a\ = 0.08, we have / ~ 6 x 10 36 
kg m 2 , so we find taml ~ 4 Gyr for Af ~ 2.6 x 10 6 s. Therefore, 
the scatter in n/Q observed for P mt > 10 days may be explained 
as a consequence of the different stellar ages. Note also that for 
WASP- 18 the high value of n/Q can be due to the very short 
orbital period which results from a very strong ti dal interaction 
in a regime with n > Q (cf. Hellie r et al ] l2009bh . 

In the case of stars of spectral types G and K we assume that 
the value of \a\ is significantly lower than in the case of F-type 
stars. This is justified because their differential rotation is lower 
than that of hotter stars, given the remarkable dependence of 
AQ on T e ff, In turn, this implies a lower toroidal field yielding 
a lower a at the same rotation period. Considering a mean value 
of |ff| ~ 0.06, we have a coronal moment of inertia I ~ 2x 10 37 
kg m 2 . For a star with the mass and the radius of the Sun, with 
Af ~ 3 x 10 6 s, this implies taml ~ 500 Myr. Such a timescale 
corresponds to that of the initial fast angular momentum loss 
occuring on the main sequence during the transit i on bet ween 
the two braking sequences introduced by iBarnesI d2003l) . i.e., 
from the so-called convective to the interface sequence. Note 
that for a mid-G-type star without a close-in planet such a tran- 
sition occurs in 100 - 300 Myr. Therefore, the effect of a hot 
Jupiter is that of slowing down the initial angular momentum 
evolution of G stars by a factor of ~ 2-5. The same is true also 
for K-type stars, but, since their transition from the convective 
to the interface sequence takes longer (« 500-800 Myr), the 
effect of the close-in planet is less important. 

In conclusion, in the case of a G- or K-type star accompa- 
nied by a hot Jupiter, we expect a significant slowing down of 
the initial phase of its rotational braking, particularly when its 
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initial rotation period is shorter than 8-10 days and the star 
has a sizeable photospheric azimuthal field component. When 
its rotation period becomes longer th an ^10 days, it s toroidal 
field component declines steeply (cf., iPetit et al]2008l) leading 
to a decrease of \a\ and a remarkable increase of the angular 
momentum loss rate. In this phase, the rate of angular momen- 
tum loss might become similar to that of stars without planets 
and the subsequent evolution could not be remarkably affected 
by the presence of a hot Jupiter, i.e., the star would continue to 
spin down according to the usual Skumanich la w charac t eristic 
of sta rs on the so-called interface sequence of iBarnesI d2003l 
2007'[ cf. Sect. |2J. Considering the different ages of planet- 
harbouring stars and their different initial rotation periods, we 
may explain the larger dispersion of n/Q observed in stars with 
T eS £5800K(cf.Fig.[D. 

In our treatment of the main-sequence spindown we have 
assumed that a star is braked as a rigid body (cf. Eq.fToli. This 
hypothesis is adequate in the present case because our braking 
time scales taml are generally longer than the time scale for an- 
gular momentum exchange between the radiative interior and 
the outer convection zone which evolutionary models of stellar 
rotation set at ^ 100 Myr on the main sequence (cf. iBouvierl 



2008; Irwin & Bouvier 2009). For the same reason, the tidal 



synchronization time should be computed by considering the 
spin-up of the whole star, as we did in Sect. |21 not just of its 
convection zone. 



4.3. A tentative comparison with observations 



In the framework of a Skumanich-type braking law, Barnes 
(2007) provides an empirical formula to estimate the age of 
a main-sequence star from its rotation period and colour in- 
dex. We apply it to HD 149026, HAT-P-1 and WASP-15 to 
test the predictions of our model for stars with a rotation pe- 
riod P mt > 10 days and 6000 < T eW < 6300 K. These three 
systems have been selected because they have a tidal synchro- 
nization time at least 3 times longer than their maximum es- 
timated ages, in order to exclude tidal effects on their angular 
momentum ev olution. Their ages, as estimated with Eq. (3) of 
Barnesl (120071) are 2.2, 2.0 and 6.0 Gyr, respectively. They are 
all within the range of ages estimated by isochrone fitting, as 
reported in Table [3] For the first two stars, the gyrochronol- 
ogy ages are close to the lower limit given by isochrone fit- 
ting, while for WASP-15, the gyro age is close to the isochrone 
upper bound. Therefore, this preliminary comparison suggests 
that some reduction of the angular momentum loss rate may 
still be induced by a close-in massive planet when P rot > 10 
days and 6000 < T e ff < 6300 K, at least in some cases, al- 
though this needs to be confirmed by a larger sample of sys- 
tems. Note that in Sect. [2] we found a similar result based on a 
greater sample of stars providing us with significant statistics. 
However, in that case we took into account the evolution of an- 
gular momentum and the tidal effects in separate analyses in 
order to have a significant sample in both cases. Now, we have 
considered the evolution of the angular momentum of stars se- 
lected to have negligible tidal effects, which severely restricts 
our sample. 



A major limitation of the present approach is that stellar 
ages derived from isochrone fitting are highly uncertain, espe- 
cially for stars with M < 1 M . Therefore, better age estimates 
are needed, such as those derived for open cluster members. 
Searches for transiting planets in open clusters have just begun 
and it is hoped that they may contribute to cla rify this issue 
(e.g., iMontalto et al.ll2007tlHartman et al.ll2009bl) . 



5. Conclusions 

We have analysed the rotation of stars harbouring transiting hot 
Jupiters and have found a general trend toward synchronization 
with increasing effective temperature. Stars with > 6500 K 
are synchronized or have a rotation period close to twice the 
orbital period of their planets (n/Q. ^ 1 or 2, respectively), 
while those with 6000 < T eff < 6500 K have n/Q a 1 or 2 only 
for P 10t < 10 days. Stars with ^ 6000 K generally show 
rotation period remarkably longer than the orbital periods of 
their planets. 

We conjecture that planet-harbouring stars are borne with 
circumstellar discs in which hot Jupiters form and migrate in- 
ward while the disc locks the rotation of the star. Depending 
on the magnetic field strength in the inner region of the disc, 
two different migration scenarios are possible, leading to a state 
with n/Q a 1 or n/Q a 2, respectively. When the discs disap- 
pear, most of the stars with M < 1.5 M are still contracting 
toward the ZAMS, so their rotation accelerates owing to the re- 
duction of their moment of inertia. Nevertheless, we conjecture 
that the synchronization between stellar rotation and planetary 
orbit is maintained throughout the final phases of the PMS evo- 
lution by the s trong coupling prov ided by a magnetocentrifugal 
stellar wind dLovelace et al.ll2008l) . Stars with M > 1.4-1.5 
M and very long lived discs (w 15 Myr) may arrive on the 
ZAMS while still locked to their discs, thus starting their main- 
sequence evolution in a status with n/Q a 2. 

Once a star has settled on the ZAMS, its rotational evolu- 
tion is ruled by the angular momentum loss from its corona. We 
assume that stars accompanied by close-in giant planets have 
a coronal magnetic field dominated by closed field lines, so 
that most of their angular momentum loss occurs through erup- 
tions similar to the solar coronal mass ejections rather than via 
a continuously streaming stellar wind. This peculiar configu- 
ration is induced by the steady motion of the planet through the 
stellar corona which reduces the magnetic helicity of the coro- 
nal fie ld leading to a predominance o f closed magnetic loops 



(e.g., lLanzal l2009; Cohe n et al. 2009). Using a simple linear 



force-free field, we estimate the angular momentum loss rate 
for different field geometries characterized by different values 
of the force-free parameter a. We find that the angular momen- 
tum loss decreases by two orders of magnitude when \a\ ranges 
from ~ 0.05 to ~ 0.2. 

If \a\ ~ 0.15 - 0.2 is characteristic of F-type stars with 
r eff ^ 6000 K and ZAMS rotation periods P mt 5 10 days, their 
rotational evolution requires timescales of the order of 30 - 
100 Gyr, that is those planetary systems would be characterized 
by an almost constant distribution of spin and orbital angular 
momentum all along their main-sequence lifetime, with their 
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present status reflecting their angular momentum distribution 
on the ZAMS. 

On the other hand, F-type stars with a rotation period ini- 
tially longer than ~ 10 days are characterized by a smaller 
value of a, say, ~ 0.08 - 0.1, leading to a greater angular 
momentum loss during coronal mass ejections. Their spin is 
expected to evolve on a timescale of ~ 4 - 7 Gyr, leading to 
some spreading in the distribution of n/Q in the effective tem- 
perature range 6000 - 6500 K as a consequence of the differ- 
ent ages of the stars. Later-type stars are characterized by still 
smaller values of a, i.e., ~ 0.05, leading to shorter braking time 
scales. Therefore, the angular momentum evolution of planet- 
harbouring G- and K-type stars should not be dramatically dif- 
ferent from that of stars without close-in massive planets (cf. 
Sect. 0. However, a reduction of the angular momentum loss 
rate by a factor of 2 - 5 may still be caused by planets around 
young, rapidly rotating (f rot $ 10 days) stars. 

Such predictions can be tested by increasing the sample of 
F, G and K stars with known hot Jupiters, especially in open 
clusters of different ages allowing us to compare the rotational 
evolution of coeval stars with and without close-in planets. 
However, since open cluster members are usually faint, this 
requires dedicated programs to be conducted with large tele- 
scopes to reach the necessary photometric and radial velocity 
precisions. 

Asteroseismology can provide stellar ages with an accuracy 
of =s 10 percent of th e total stellar main-sequence lifetime (e.g., 
Kjeldsen et all 2009). but the internal chemical composition of 
planet-hosting stars may differ from t hat of their surface laye rs 
inducing systematic errors (cf.,e.g., Baz ot & Vauc"laull2004h . 

In principle, spectropolarimetric techniques can be applied 
to derive the value of the parameter a in stars harbouring 
hot Jupiters, provided that they rotate sufficiently fast (cf. 
Sect. 13.11 ). This should allow us to test our theory in the case 
of individual objects, at least those with a sufficiently rapid ro- 
tation. 

The possible effect of hot Jupiters on stellar angular mo- 
mentum l oss must b e taken into account when interpreting the 
results of Pont (2009) in the sense that they could not necessar- 
ily provide evidence that tides are ruling the spin evolution in 
stars with close-in planets. It is more likely that both tides and 
the effects discussed in this paper are simultaneously at work to 
affect the distribution of angular momentum and its evolution 
in stars harbouring hot Jupiters. 

Finally, we note that gyrochronology may not be suitable to 
estimate the age of late-type stars with close-in giant planets, 
especially if they have T e ff > 6000 K and/or are rotating with 
a period shorter than ~ 10 days, because their rotational evolu- 
tion can be remarkably different from that of stars without hot 
Jupiters. 
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